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THERMAL  CONDUCTIVITY  OF  SUPERCONDUCTORS 
INTRODUCTION 

The  discovery  of  electrical  superconductivity  In  1911  triggered 
experiments  designed  to  determine  whether  the  phenomenon  occurs  in  heat 
transport.  The  first  rf,sults^  showed  that  instead  of  thermal  super¬ 
conductivity,  the  heat  conduction  in  the  superconducting  state  Was  lower 
than  in  the  normal  state.  Since,  in  metals,  in  the  normal  state,  -the 
conductivity  is  composed  of  the  electronic  component  plus  the  lattice 
component,  it  may  be  written  as 


In  the  superconductive  state,  we  may  write 


K,  =  ,  +  K  ,,  ,  . 

s  e/,s  ph,s 


(2) 


Here  the  subscript  n  stands  for  the  normal  state,  s  for  the  superconducting 
state,  e/  for  the  electronic  contribution  to  the  thermal  conductivity,  and  ph 
the  lattice  contribution.  Since  the  "suoerelectrons"  which  carry  the  re- 
sistanoeless  current  in  a  superconductor  cannot  be  expected  to  carry  the 
thermal  current,  K  .  would  be  expected  to  be  smaller  than  k  .  J''  Likewise, 
since  the  supereiectfons  move  without  friction  against  the  lattice,  they  pre¬ 
sumably  do  not  contribute  to  the  scattering  of  the  lattice  waves,  hence  >'  u  5 
might  be  expected  to  be  larger  than  k  ^  Therefore,  since  the  heat  is^  ' 
carried  only  by  that  fraction  of  the  electrons  which  are  still  normal  and  by 
the  lattice  component,  would  be  expected  to  be  smaller  than  k  •  because 
the  lattice  contribution  is  small  in  both  states.  Qualitatively,  this  ex¬ 
plains  the  results. 

A  complete  understanding  of  the  thermal  conductivity  in  superconductors 
can  only  be  achieved  as  a  consequence  of  a  detailed  microscopic  theory  of 
superconductivity.  The  thermodynamic  and  electrodynamic  behavior  of  super¬ 
conductors  has  led  to  a  working  hypothesis  which  has  been  remarkably  success¬ 
ful  in  a  rough  interpretation  of  the  observed  effects,  although,  in  Its  crude 
form.  It  cannot  have  any  physical  significance.  This  is  the  two-fluid  model 
In  which  the  electron  fluid  is  regarded  as  a  completely  Interpenetrating  mix¬ 
ture  of  a  normal  and  a  superconductive  constituent. 

The  Two-Fluid  Model 

The  two-fluid  model  does  not  explain  the  phenomena  of  superconductivity 
but  is  a  convenient  scheme  for  their  description.  It  assumes'*  that  a  frac¬ 
tion  (1-X)  of  ti  3  Fermi  surface  is  modified;  the  electrons  on  the  surface 
condense  into  a  lower  state.  It  also  assumes  that  the  electrons  in  those 
modified  regions  cannot  be  thermally  excited,  though  the  fraction  of  the  Fermi 
surface  thus  affected  is  a  function  of  temperature  and  increases  with  decreas¬ 
ing  temperature.  .The  superconductive  regions  can  be  oriented  so  as  to  yield 
a  supercurrent,  it  has  been  shown  experimentally  th&t  the  decrease  in  total 
entropy  is  a  result  of  the  growth  of  the  superconductive  concentration  and 
that.  In  fact,  the  entropy  of  this  constituent  is  zero  at  all  temperatures 


below  the  critical  temperature  T^.  Since  the  superconductive  regions  con¬ 
tribute  zero  to  the  entropy,  there  being  no  thermal  excitation,  the  Thomson 
coefficient  for  supercurrents  Is  zero,  and  the  entropy  Is  less  than  It  would 
be  In  the  normal  state.  Similarly,  the  electronic  thermal  conductivity  Is 
reduced,  since  only  the  normal  fraction  of  the  Fermi  surface  contributes 
toward  It.  Such  conclusions  led  to  the  postulate  that  there  Is  an  energy  gap 
In  the  electron  spectrum  of  the  metal”  which  Is  roughly  coincident  with  the 
Fermi  energy.  The  Bardeen,  Cooper,  Schrleffer  (BCS)  theory”  and  some  experi¬ 
mental  data^  seem  to  favor  such  a  model.  On  the  other  hand,  the  lattice 
component  of  thermal  conductivity  Is  enhanced,  since  the  lattice  waves  can 
only  be  scattered  by  the  electrons  on  the  normal  fraction  of  the  Fermi  surface. 
Another  explanation  Is  that  In  the  superconducting  state  the  lattice  thermal 
conductivity  Is  Increased  relative  to  that  In  the  normal  state  because  the 
anergy  gap  In  the  elactronlc  spectrum  leads  to  an  Increase  In  the  relaxation 
time  for  phonons.  The  energy  gap  confers  on  the  metal  an  aspect  which  Is  not 
too  different  from  that  of  a  dielectric  crystal.”  Owing  to  the  small  size  of 
the  gap,  this  behavior  can  not  be  seen  except  at  very  low  temperatures.  This 
feature  of  superconductivity  (dieIectrIc-1 Ike  behavior)  becomes  very  convinc¬ 
ing  In  the  heat  conduction.  However,  because  of  the  smallness  of  the  lattice 
term,  the  electronic  contribution  (though  rapidly  falling)  will  dominate  the 
thermal  conductivity  Just  below  T  . 

C 

Assume^*”  that  the  modification  of  the  superconducting  fraction  of  the 
Fermi  surface  has  associated  with  it  a  latent  heat;  thus  the  Helmboltz  free 
energy,  usually  of  the  form 

F--I/2yT”,  (3) 

Is  now  modified  to 

F  -  -  1/2  ya  yT^  -  (l  -  y)R  .  (k) 

The  first  term  Is  the  contribution  from  the  normal  region,  the  second  term 
Is  the  contribution  from  the  superconducting  region,  y  and  0  are  constants 
characteristic  of  their  respective  regions.  The  latter  term  Is  due  to  the 
latent  heat  since  there  Is  no  continuum  of  states  available  for  thermal  ex¬ 
citation,  The  condition 

(aF/ay)^  -  0  (5) 

gives  X  ””  ”  Function  of  T.  It  Is  not  possible  to  assume  simply  a  ■  I, 
as  It  would  be  In  the  absence  of  Interaction  between  the  n-  and  s-  regions, 
since  this  would  not  generally  satlsfy(5)« If  one  takes  a  ■  1/2,  the  observed 
thermodynamic  properties  ere  approximately  reproduced,  I.e. 

y  -  (yT*/1*R)*  -  (T/T^)\  (6) 

where  T  Is  the  transition  temperature,  for  which  y-  I,  The  specific 

C 


2 


heat  per  unit  volume  Is 


C  =  -  T  (d^F/dT^)  =  lyd^/Tj') .  (7) 

but  the  specific  heat  due  to  the  thermal  excitation  of  the  electrons  in  the 
normal  region  is 

C  -  -  T(d^F/dT=)  -  y(TVT  •■’).  (8) 

n  X  c 

The  difference  C  -  C  is  ascribed  to  the  change  of  energy  as  electrons 
change  their  phase,  f.e.  from  being  In  an  s-reglon  to  being  In  an  n-reglon. 

The  thermal  conductivity  of  normal  metals  can  be  written 
K  =  1/3  C  V  "  (9) 

where 

C  =  ( 1 0) 

is  the  electronic  specific  heat,  v  is  the  velocity  of  electrons  of  Fermi 
energy  (assumed  isotropic)  and 

/  =  VT  ( 1 1 ) 

is  the  effective  electron  mean-free-path.  The  two-fluid  model  has  been 
applied^®  to  the  thermal  conductivity;  the  electronic  thermal  conductivity 
in  the  superconducting  state  differs  from  that  in  the  normal  state  because  C 
and  possibly  t  are  altered.  For  C  we  should  now  use  C  ,  because  we  are  con¬ 
cerned  with  the  transport  of  energy  by  electrons  which  Pemain  normal  when 
passing  along  a  temperature  gradient,  rather  than  with  a  change  of  energy 
due  to  a  change  of  phase.  Thus, 


’*e'',s^’^e/, 


=  (T/T.)^('  /! 


(12) 


It  should  be  noted  that  if  we  had  used  the  C  defined  by  (7)  instead  of  C 
defined  by  (8)  to  calculate  h  ,  ,  there  would  have  been  a  discontinuity^ 

at  the  transition  tempera ture^lrt’the  curve  of  h  .  vs  T,  similar  to  that  ob¬ 
served  for  the  specific  heat,  because  /  (T  )  mult  equal  f  (T  ) , 

s'  c  ^  n  c 

The  behavior  of  *  requires  some  assumption;  '  could  differ  from  i 
because,  when  an  election  is  scattered  from  a  state^in  the  n-region,  thePe 
are  fewer  final  states  available  to  it,  since  the  states  in  the  s-region 
are  modified.  The  ratio  '  /•'  should  depend  on  the  mechanism  of  scattering. 
It  was  suggested^®  that  fo?  spattering  by  static  imperfections 


s 


2  »,„/(!  ♦  k) 


(13) 


7 


so  that  i  >  This  form  would  be  appropriate  If  scattering  were  Isotropic 
and  If  a  fraction  x/2  of  the  possible  final  states  were  blocked  by  electrons 
In  the  s-states.  Thus,  If  the  thermal  resistance  Is  mainly  due  to  static  im¬ 
perfections,  the  fractional  change  of  the  conductivity,  on  passing  from  the 
normal  to  the  superconducting  state  Is 


“o  n^«  ,  "  •‘o  n  “  2(T/T  )VI  +  (T/T  )  ^ 

OpSOp^  c  c 


Consider  now  the  case  where  the  thermal  resistance  Is  mainly  due  to 
lattice  waves:  at  low  temperatures  an  electron.  Interacting  with  a  phonon, 
does  not  change  Its  "horizontal"  position  on  the  Fermi  surface  (from  a  po¬ 
sition  or  region  of  large  concentration  to  that  of  a  lower  one)  by  a  large 
amount  (see  reference  II  p.9  and  Fig.  1).  it  thus  appears  that  an  electron 
In  the  n-reglon  will.  In  the  majority  of  cases,  remain  In  the  n-region  after 
an  Interaction,  so  that  i  •I  ,  Hence  the  change  In  the  ideal  (intrinsic) 
thermal  resistance  ^ 


It  would  hardly  be  expected  that  (l4)  and  (15)  would  give  the  temperature 
dependence  of  the  ratios  f  and  g  exactly,  although  one  would  expect  these 
equations  to  give  at  least  a  qualitative  description  of  their  variation.  We 
shall  see  later  that  this  Is  so  for  (|4)  but  not  for  05). 


The  two-fluid  model,  In  the  explicit  form  given  above,  reproduces  the 
thermodynamic  properties  only  at  temperatures  above  about  T  /2;  at  lower 
temperatures,  the  specific  heat  decreases  exponentially  wItR  decreasing  tem¬ 
perature,  and  In  view  of  (9)  one  would  expect  k  ,  to  behave  similarly.  It 
should  be  noted  that  at  T  ,  the  superconductive  nftSt  conductivity  k  breaks 
away  from  h  .  The  change*'may  be  sudden  or  gradual^  Fig.  I,  but  no^discon- 
tlnulty  In  h  Is  observed  at  T  even  though  this  effect  has  been  carefully 
investigated.^  However,  at  T^  there  Is  a  discontinuity  In  the  specific  heat. 
Fig.  2,’®  This  jump  Is  of  magnitude 


(16) 


[compare  with  equation  (lO)].  Measurement  of  y  provides  direct  information 
about  N(E_) ,  the  density  of  states  at  the  Fermi  level.  At  lower  temperatures 
the  energy  gap  ^(T)  tends  to  dominate  the  specific  heat,  but  a  simple  formula 


I  Ike 


Ce^  s  ^  exp(-MT)/T)  (17) 

Is  not  adequate,  until  t  reaches  Its  limiting  value  which  is  expressed 

A  ea  1 . 76  J<  T  (18) 

o  —  c 

where  A  Is  the  value  of  the  energy  gap  at  T  -  0  and  J<  Is  the  Boltzmann 
constant, 
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The  Ratio  y.  /h 
s  n 


Since  the  electrons  In  the  superconductive  region  of  the  Fermi  surface 
cannot  be  thermally  excited  Into  a  continuum  of  states,  it  follows  that 
lattice  waves  can  be  scattered  only  by  the  electrons  of  the  normal  region, 
and  not  by  those  of  the  superconductive  regions  of  the  Fermi  surface.  Thus, 


(T/T^) 


(19) 


where  W_  is  the  thermal  resistance  due  to  the  conduction  electrons,  and  if 

K  .  Is  limited  by  the  interaction  with  conduction  electrons 
ph 


K  .  ./k  . 
ph,s  ph,n 


1/x 


(20) 


However,  the  lattice  resistance  due  to  phonon-phonon  interactions  should  be 
unchanged  by  the  transition  from  the  normal  to  the  superconducting  state.  It 
can  thus  be  seen  that  k  ^  ^  **  h  s  ^  h  n'  that,  depending  on 

the  circumstances,  mSy’Se  ei?h4?  smalle?  6r  lar§ef'^than 

The  ratio  of  electronic  thermal  conductivities  in  the  superconducting  and 
normal  states  h  ^  />;  has  been  difficult  for  superconductive  theory  to 

explain.  When  frri^urlfy’Sfcatteri  ng  is  dominant,  the  BCS  theory''  predicts^’' 
that  this  ratio  should  be  a  universal  function  of  T/T  ,  Independent  of  the 
particular  element  measured.  The  expression  is  complfcated  but  It  Is  a  very 
similar  function  to,  and  its  values  are  close  to,  those  of  the  Helsenberg- 
Koppe  formula^® ‘ 


H..,  „  •  2(T/T  +  (T/T ref  equation  (14)1.  (21) 

ei,s  eji,,n  c  c 

Both  expressions  agree  with  experiment,  and  they  have  zero  slope  at  T  »  T  . 

The  temperature  dependence  of  y  /:•  is  shown  in  Fig,  3.  However,  a  radically 
different  behavior  is  observed  in  metals  like  lead,  mercury,  and  indium, 
where  T  occurs  in  the  phonon  scattering  region  of  Fig,  4a, In  these 
elements  previous  workers  have  found  that  k  ,  /k  ,  falls  from  T  as  rapidly 
as  (T/T  ),^''  The  BCS  prediction  for  this  rSgfSn  ?S*Sompletely  at  Parlance 
with  thS  observations.^'"’ 

As  can  be  seen  from  Fig.  3,  k  should  become  extremely  small  at  suf¬ 
ficiently  low  temperatures.  However,  only  the  electronic  component  of  thermal 
conductivity  has  been  considered  in  these  calculations. This  Is  by  far  the 
dominant  mechanism  for  a  pure  metal  in  the  normal  state.  But,  it  should  be 
remembered  that  this  Is  only  the  case  because  phonon  conduction  In  metals  is 
inhibited  owing  to  the  extremely  effective  scattering  of  phonons  by  the  free 
electrons.’^,*’’’,-'’  The  removal,  on  cooling,  of  a  progressively  Increasing 
fraction  of  the  electron  fluid  from  the  thermal  distribution  in  the  super¬ 
conductive  state  not  only  decreases  the  heat  conduction  by  electrons  but  also 
decreases  the  scattering  of  phonons  by  them.  Accordingly,  phonon  conduction 
in  the  metal  will  become  the  dominant  process  at  sufficiently  low  temperatures. 
Its  behavior,  then,  will  be  closely  analogous  to  that  of  a  dielectric  crystal. 
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Thus,  at  some  temperature  below  T  ,  h  /k  must  become  larger  than  the 
function  of  Fig,  3,  So  far,  no  theoretical  attempt  has  been  made  to  determine 
this  temperature;  however,  experimental  data  suggest  that  it  will  be  below 
0.4  T  in  the  case  of  a  pure  metal.  At  sufficiently  low  temperatures,  phonon 
mean>free-paths  are  limited  by  the  size  of  the  crystal^^’^*’  and  by  its  surface 
roughness,*®**^  resulting  in  a  mean- free- path  Independent  of  temperature.  Al¬ 
so,  below  the  Debye  temperature,  6.,  phonon  velocities  are  roughly  constant 
and  the  lattice  specific  heat  varies  as  T^.  Thus,  the  thermal  conductivity 
should  become  proportional  to  T^,  In  the  case  of  good  single  crystals,  the 
factor  of  proportionality  depending  on  the  dlwneter^^  of  the  specimen  since 
the  only  relevant  process  will  be  scatter  of  phonons  on  the  walls  of  the 
specimen.  Thus  a  size  effect  in  the  heat  conductivity  of  a  superconductor 
is  expected.  Moreover,  x  ,  which  is  chiefly  limited  by  scatter  of  electrons 
on  point- imperfections  at^higher  values  of  T/T  must  become  less  sensitive  to 
these  at  lower  temperatures,  since  their  scattering  cross  section  is  small 
for  phonons.  Instead,  h  must  be  strongly  influenced  by  extended  lattice 
faults.  * 

The  relative  importance  of  conduction  by  phonons  will  be  shifted  to  high¬ 
er  reduced  temperatures  as  the  impurity  content  of  the  speciirxin  rise:.  For 
very  impure  specimens  which  have  a  low  electronic  conduction  in  the  normal 
state,  one  can  expect  values  of  h  which  approach  or  even  exceed  those  of  h  . 
The  explanation  of  the  inversion  Sf  the  k  /k  ratio  in  the  case  of  alloys 
is  thus  provided  by  a  combination  of  deprissed  electron  and  enhanced  phonon 
conduction.  The  correctness  of  this  model  is  borne  out,  at  least  qualita¬ 
tively,  by  experiments.*^  Even  so,  it  must  be  regarded  as  surprising  that, 
in  some  cases,  for  example  in  a  lead  alloy  with  10%  bismuth,  a  rapid  rise  of 

K.  occurs  already  at  T  , 
s  c 

The  opposite  case  is  given  by  a  pure  metal  with  an  intrinsically  high 
electrical  conductivity,  such  as  aluminum,  in  which  the  residual  heat  con¬ 
duction  due  to  the  normal  electrons  will  remain  the  dominant  factor  in  h 
even  at  the  lowest  temperatures  in  spite  of  an  enhanced  phonon  conductivity. 

In  such  cases  one  should  expect  a  close  adherence  to  the  semi-empirical 
formula  (f- function) ,  particularly  when,  as  in  aluminum,  T  is  well  below  the 
maximum  in  the  thermal  conductivity.  ^ 


The  Electronic  Thermal  Conductivity  of  Superconductors  with 
Strong  Electron- Phonon  Coupling 


Experiments  on  the  strong  elemental  superconductors  lead  and  mercury 
indicate  that  the  electronic  thermal  conductivity  of  these  materials  differs 
markedly  from  that  of  typical  weak  superconductors  like  tin  or  indium.  For 
lead  and  mercury®*^'’’  the  ratio,  h  /k  ,  of  the  thermal  conductivity  in  the 
superconducting  and  normal  states  'vh8n  plotted  against  the  reduced  temperature 
tA  shows  a  steep  positive  slope  of  about  5  near  T  ■  T  ,  For  tin®*  and 
indium,*^  on  the  other  hand,  the  experiments  yield  a  smaller  slope  of  about 


A  theory  of  the  electronic  thermal  conductivity  of  superconductors,  based 
on  the  quasiparticie  approximation  and  the  Boltzmann  equation  approach  of 
Bardeen,  Rlckayzen,  and  Tewordt,^^  has  been  carried  through  previously.®" 
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The  results  obtained  are  in  substantial  agreement  with  the  data  on  tin  and 
indium.  The  lack  of  agreement  of  this  theory  and  the  data  on  lead  is  not 
surprising  since,  as  shown  In  the  work  of  Schrleffer  et  al,*’’  the  strong 
electron-phonon  coupling  causes  the  quasiparticle  picture  to  be  quite  meaning¬ 
less  over  much  of  the  energy  spectrum. 

It  therefore  seemed  reasonable  to  apply  Schrleffer's  theory,  which  was 
used  in  explaining  tunneling  characteristics,  to  discuss  the  electronic 
thermal  conductivity  of  strong  superconductors  without  recourse  to  the  quasi- 
particie  approximation.*"  The  starting  point  for  this  theory  is  the  Kubo 
formula  in  which  the  thermal  conductivity  Is  expressed  In  terms  of  the  cor¬ 
relation  function  of  two  heat  current  operators.*®  This  Is  first  examined  in 
the  Hartree-Fock  approximation  in  the  Nambu***  space.  It  Is  shown  that  in  the 
Eliashberg*^  approximation  of  neglecting  the  momentum  dependence  of  the  elec¬ 
tronic  self-energy,  the  calculation  of  the  thermal  conductivity  is  reduced  to 
a  quadrature,  involving  however  the  co‘nplex  energy  gap  and  renormalization 
functions  which  are  solutions  of  the  Eliashberg  gap  equations  at  finite  tertv- 
peratures.  These  expressions  are  too  complicated  to  be  reproduced  here  but 
the  reader  is  referred  to  the  original  paper  (reference  28)  for  more  details. 
The  problem  was  also  considered  in  the  ladder  approximation  in  the  Nambu 
space.  A  generalized  Boltzmann  equation  was  derived  which  Includes  corrections 
to  the  Hartree-Fock  approximation  corresponding  to  the  replacement  of  the  scat¬ 
tering  by  the  transport  cross-section.  This  treatment  shows  that  the  standard 
Boltzmann  equation  for  superconducting  quasi-particies  is  obtained  In  the  weak- 
coupl ing  limit. 

The  thermal  conductivity  of  pure  superconducting  lead  has  bean  considered 
to  be  anomalous  for  many  years.''*  The  experimental  results  are  summarized  in 
Fig.  5.  As  was  mentioned  above,  the  curve  of  )-  /h  vs  T/T  for  typical  weak- 
coupling  superconductors  tin  and  indium  has  a  smalY  limiting  slope  of  about 
1.3.  in  the  case  of  the  strong-coupling  superconductors,  lead  and  mercury, 
the  decrease  in  the  thermal  conductivity  is  steeper,  i.e.  the  limiting  slope 
for  lead  is  about  9.®'^ 

The  strong-coupling  superconductors  are  characterized  by  large  electron- 
phonon  matrix  elements,  and  by  peaks  at  low  energies  in  the  density  of  phonon 
states  in  which  the  electrons  are  coupled.  These  distinguishing  character¬ 
istics  have  unambiguously  been  shown  to  be  responsible  for  the  anomalously 
large  values  in  lead  and  mercury  of  the  ratio  of  the  energy  gap  at  0  °K  to  the 
critical  temperature,  and  for  the  anomalous  thermodynamic  properties  of  these 
metals.*  It  has  often  been  speculated  that  the  smaller  thermal  conductivity 
of  these  s trong- coup  1 ing  superconductors  is  another  consequence  of  their  un¬ 
usual  electron-phonon  interactions.  It  was  not  clear,  heretofore,  how  this 
explains  the  great  reduction  in  thermal  conductivity. 

The  calculation  is  based  on  the  foregoing  theory,*"  and  the  general 
theory  is  supplemented  by  a  specific  model’''  for  the  phonon  spectrum  and  the 
electron-phonon  coupling  constants  in  lead.  The  phonon  spectrum  used  is 
shown  in  Fig,  6.  it  is  found  that  even  near  the  critical  temperature,  long- 
lived,  parti cle- 1 i ke  excitations  exist  for  the  energies  important  in  thermal 
conduction,  in  this  quasiparticle  limit,  the  general  formula  obtained  in 
reference  28  (equation  2,17)  reduces  to  (neglecting  the  effect  of  scattering- 
in  terms) 
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d(j)  mruf  A|  T)  ]^/Z,  r^(u),T)  , 


(22) 


« 


K 


S 


A/r  r 

a,  (A,,T) 


where 


A  «  N(0)  k 


(23) 


with  N(0)  the  density  of  states  at  the  Fermi  surface  for  one  spin,  Vp 
the  Fermi  velocity,  uo  the  angular  frequency  of  the  phonons.  A,  ((d,T) 
the  real  part  of  the  Eiiashberg  gap  function®® Z  the  renormalization 
function,  r  (uj,T)  the  quasiparticle  1 i fetime  which  is  related  to  the 
parameters  of  the  Eiiashberg  theory  according  to®"**'*® 

u,Z,r(ti))  -  22a (uis  -  A,"")  -  2A,AaZ,  ,  (24) 

J<  is  Boltzmann's  constant,  and  T  the  absolute  temperature.  Equation  (22) 
has  the  same  general  form  as  is  obtained  from  a  phenomenological  Boltzmann 
equation,^®*®*®  But,  there  is  an  Important  difference  because  the  virtual 
effects  of  phonons  in  causing  the  superconducting  transition  have  not  been 
approximated  by  a  model  potential,  but  have  been  treated  in  the  same  way 
as  the  real  transitions  which  scatter  quasiparticles.  As  a  result,  the  large 
value  of  the  energy  gap  (in  units  of  k  T  ) ,  is  in  principle  and  practice 
contained  In  (22).  ^ 

In  order  to  bring  out  the  physical  origins  of  the  large  limiting  slope 
of  the  reduced  thermal  conductivity  of  lead  (h  /k  vs  T/T  ),  it  is  helpful 
to  have  a  formal  expression  for  this  slope.  Such'^an  expression  follows  at 
once  from  (22),  The  thermal  conductivity  in  the  normal  state  Is  obtained 
from  this  equation  by  setting  A  equal  to  zero.  Forming  the  ratio  x  /x  and 
taking  the  derivative  at  the  critical  temperature,  one  obtains  * 

d/dt(H5/x^)  -  -  l/23^®(Su,''/3t)' 
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o 

sech®  1/23  >i}  (25) 

where  3  ■  ^/J<T,  t  “  F  2,  r,  and  the  temperature  derivative  of  A,®  is 

taken  outside  the  integration  because  it  is  essentially  constant  in  the 
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relevant  region  of  lu.  The  following  three  factors  appear  to  be  responsible 
for  the  large  slope  In  lead  as  contrasted  with  weak-coupl ing  materials: 

1,  Tne  larger  value  of  the  ratio  2^»(0)/J<  T  (4,3  for  lead  as  opposed  to  3.5 
for  materials  well  described  by  the  DCS  theSry)  has  as  its  corollary  a  larger 
value  of  the  slope  -  ")  (1^.1  for  lead,  9.4  for  materials  follow¬ 

ing  BCS) .  The  more  rapid  openinn-up  of  the  energy  gap  in  lead  means  physical¬ 
ly  that  the  heat  carryin'  gtias ipar! iclcs  are  more  rapidly  frozen  out.  This  is 
the  most  obvious  cause  o'"  ti'e  reduced  thermal  conduct  i  i  ty,  but,  taken  by  it¬ 
self,  it  does  not  suffice  to  explain  the  large  effect,  2,  The  quasiparticle 
lifetime  decreases  more  rapidly  with  frequency  in  lead  than  in  weak-coup- 
ling  materials.  This  e.ffect  has  Its  origin  in  the  small  density  of  low-fre¬ 
quency  phonons  in  lead.  The  ratio  of  integrals  multiplying  ^/''t(A|)”'  In  (25) 
is  the  larger  the  more  rapidly  decreases  with  frequency.  In  loose  physi¬ 
cal  terms  one  can  say  that  in  all  materials,  the  advent  of  the  energy  gap 
suppresses  the  carriers  that  are  most  weakly  damped,  and  are  thus  most  ef¬ 
ficient  in  carrying  energy.  This  theory  indicates  that  in  lead  this  supres- 
sion  is  particularly  effective.  For  the  ratio  of  integrals  mentioned,  a  value 
of  about  1,1  is  found.  For  the  model  of  Debye  phonons  and  "jellium"  matrix 
elements  worked  out  by  Tewordt,^"*’  the  ratio  is  aljout  0.5.  3.  The  ratio 

”  (u:,T)/  '■  ((u,T)  decreases  for  lead  v^hen  T  decreases  below  T  so  that  the 
last  ''atio  of  integrals  in  (25)  is  positive.  The  sign  of  thfs  term  appears 
to  be  connected  with  the  coherence  factors  v/hich  go  into  a  calculation  of  the 
rela”ation  rate  for  a  auasipartide  in  a  superconductor.  In  this  model,  the 
dominant  relaxation  process  is  one  in  which  two  quasiparticle  excitations 
annihilate,  emitting  a  phonon.  This  is  labeled  process  3  In  Fig,  7,  This 
gives  a  positive  sign.  The  other  two  kinds  o'  processes,  labeled  1  and  2 
in  Fig,  7,  scattering  of  quasiparticles  witli  phonon  emission  and  absorption, 
give  the  opposite  sign.  The  value  for  this  term  ''the  last  ratio  of  integrals 
in  (25)1  is  calculated  to  be  3.5.  V/orklng  backwards  from  the  final  slope 
obtained  in  reference  26b,  it  can  be  concluded  that  for  the  model  used  in  this 
reference,  the  ratio  is  negative  and  approximately  -0,9, 

Although  no  one  of  the  three  factors  discussed  is  large  enough  to  account 
for  the  effect,  taken  together  tioy  change  th.e  slope  1,6  obtained  in  Te- 
wordt's^''’°  model  to  the  large  'aliie  II.  In  Towordt's  model  the  virtual  pro¬ 
cesses  are  accounted  for  by  the  BCS  rx)del  and  the  real  processes  by  a  Debye 
spectrum  of  longitudinal  phonons  coupled  to  the  electrons  by  "iellium"  matrix 
elements. 

The  model,  although  containing  the  effects  discussed  above,  has  certain 
weaknesses.  The  absolute  value  of  the  thermal  conductivitv  at  the  critical 
temperature  has  approximately  the  correct  value.  At  lower  temperatures,  how¬ 
ever,  the  complete  absence  of  lowTrequency  phonons  will  result  in  the  thermal 
conductivity  in  the  normal  state  rot  approaching  the  T"^  increase  of  the  Bloch 
theory,  but  instead,  increasing  exponentially.  In  spite  of  Such  weaknesses, 
the  basic  conclusions  of  this  thcorv  appear  sound. 

This  theory  is  based  on  a  pre  ious  theory"*'  whict)  was  motivated  by  a 

Note  in  Fig.  6  that  there  are  no  longitudinal  phonons  below  7  meV  and  no 

phonons  whatever  below  2.15  meV(mi 11 i-electron  volts). 


feeling  that  the  quasiparticle  approximation,  in  the  sense  of  iifetimes  being 
small  compared  to  excitation  energies,  might  break  down  for  thermai  conduc¬ 
tivity  in  lead.  The  present  theory  indicates  that  no  such  breakdown  occurs. 
However,  the  virtual  effects  of  high  energy  phonons  are  Important  for  thermal 
conductivity.  These  are  consistently  treated  by  the  previous  theory*®  and  not 
by  a  phenomenological  Boltzmann  equation.  The  results  are  in  fair  agreement 
with  experiment. 

The  difference  in  behavior  between  the  strong-coupling  and  the  weak- 
coupling  superconductors  is  clearly  apparent  from  Fig,  I;  the  curve  depart¬ 
ing  from  H  abruptly  in  the  case  of  lead  and  gradually  in  the  case  of  tin. 

It  can  a1s8  be  seen  that  the  transition  temperature  T  ,  for  lead.  Is  above  the 
maximum  In  x  and  for  tin,  it  is  below  it.  This  means  that  the  scatter  of  the 
normal  elect?ons,  at  the  onset  of  superconductivity,  is  mainly  by  phonons  In 
lead  and  by  impurities  in  tin.  In  both  cases,  however,  is  lower  than 
for  the  whole  range  of  superconduct  ivlty. 

Thermal  Conductivity  in  the  Superconductive  State 


The  usual  procedure  Is  to  measure  the  thermai  conductivity  of  super¬ 
conductors  both  In  the  normal  and  the  superconductive  states.  The  former 
measurements  are  made  in  a  magnetic  field  larger  than  that  needed  to  quench 
superconductivity,  and,  if  necessary,  are  reduced  to  zero  field  strength  by 
extrapolation.® 

Observations  on  the  thermal  conductivity  in  the  superconducting  state 
can  be  classified^  into  (a)  cases  where  h  .  (  h  .  and  W.  )  W  ,  (b)  cases  where 
H  t  <  K  .  end  W  )  W . ,  (c)  cases  where  h  n?§Iigibli  in  ?he  normal  state 

bBl  appreciable^! n  the  superconducting  s*2ate,  and  (d)  cases  where  x  .  is  ap¬ 
preciable  both  in  the  normal  and  In  the  superconducting  states,  ThBre  are, 
of  course,  cases  intermediate  between  any  of  the  above  classes,  and  their  in¬ 
terpretation  is  correspondingly  uncertain.  W  and  W.  are  the  residual  and 
Ideal  (intrinsic)  thermal  resistivities,  respicti vely.  Thus,  as  In  the  case 
of  thermal  conductivity  in  normal  metals,  the  thermal  conductivity  of  super¬ 
conductors  can  be  discussed  in  terms  of  the  scattering  of  electrons  by  lattice 
waves,  i.e.  the  ideal  thermal  resistance  W. ,  the  elastic  scattering  of  elec¬ 
trons  by  impurities,  i.e.  the  residual  thermal  resistance,  W  ,  and  the  lattice 
component  of  thermai  conductivity.  ° 


The  ideal  Resistance  in  the  Superconducting  State 


To  realize  the  condition  W  )  W  below  T  ,  the  specimen  must  be  very 
pure  and  the  transition  temperature  Reasonably  high.  This  condition  has,  so 
far,  been  fulfilled  only  in  lead  and  mercury.  The  thermai  conductivity  of 
lead  is  shown  in  Fig,  lb;  Figs,  5  and  8  show  h  /h 
mercury,  ®  " 


vs  T/T^  for  both  lead  and 


From  the  reasoning  used  to  develop  (15),  H  should  be  independent  of  T 
just  below  T  ,  until  scattering  by  Imperfections  becomes  important.*  in  gen¬ 
eral,  H  should  be  of  the  form 


1/k  «  W  -  W,  +  W 
s  s  l,s  o,s 


(26) 
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where 


W.  =  W  (T  ) 

I  ,s  I ,n  c 

is  independent  of  temperature  and  ^  is  related  to  ^  by  (\k),  white 

w  « .  (28) 

o,n 

For  this  discussion,  it  is  I>riportant  to  note  that  W  ^  should  increase 
monotonical ly  with  decreasing  temperature,  and  so  sRduld  W^, 

The  observed  behavior  of  h  does  not  conform  to  these  predictions. 
Immediately  below  T^, 

«  T"' ,  (29) 


so  that 


g  -  (T/T^)^  (30) 

in  contrast  to  (15).  At  lower  temperatures  does  not  decrease  steadily 

v;i  th  decreasing  temperature,  but  increases  again  and  then  decreases  at  a 

temperature  such  that  W  is  comparable  to  W.  , * 

o,n  I  ,n 

There  are  at  least  three  interpretations  of  the  observed  behavior  of 
K  The  one  considered  most  likely  by  Kiemens'*  is  that  W.  ^  is  approxl- 
mitely  described  by  (15)*  but  for  some  unknown  reasons  whiih  are  presumably 
outside  the  scope  of  the  two-fluid  treatment,  (15)  does  not  describe  g  im¬ 
mediately  below  T  ,  Just  below  T  ,  the  actual  g-function  is  smaller  than 
(15),  but  not  by  I  very  large  amount  '"the  observed  values  of  h  are  not  less 
than  3/^  k  (T  )  at  the  minimum"';  at  lower  temperatures  k  Increases  again  to 
its  theoreticSl  value  i-  (T  ) ,  then,  at  still  lower  temperatures,  decreases 

monotoni cal ly,  due  to  W  , 

0,5 

Then  there  is  the  interpretation  that  (30)  holds  over  a  wide  range  of 

temperatures,  but  that  the  ideal  and  imperfection  resistances  do  not  combine 

additively  as  in  (26),  so  that  W  is  not  a  monotonic  function,  even  though 

W  and  W.  are.  This  viewpoint  raises  two  questions:  (1)  Why  should 

be  so'iTiuch  smaller  than  \  ?  and  (2)  Why  should  (26)  break  down  so 

vfdfently  at  intermediate  tempefatures  that  W  decreases  with  decreasing 

temperatures,  even  though  both  W.  and  W  increase? 

t ,s  o,s 

The  third  interpretation  ascribes  the  maximum  in  k  to  an  enhanced 
lattice  thermal  conductivity  in  the  superconductive  state,  but  the  difficulty 
is  that  a  rather  large  amount  of  lattice  conductivity  is  needed,  i.e.  k 
{h  °K)  r-t  1  watt/cm-deg.’ By  direct  evaluation  of  a  theoretical  formu-** 
la''^,  an  estimate  for  lead  was  obtained  which  implies  that  ;•  (U  °K)  <  10“® 
watt/cm-deg.  Also,  measurements''''  on  lead  showed  that  h  wa?"llmlted  only  by 
boundary  scattering  at  the  lowest  temperatures.  An  upper  limit  of,  say  5%, 
can  be  put  on  the  resistance  of  all  other  phonon  mechanisms  at  1  °K,  If  this 
57  were  all  due  to  the  scattering  of  phonons  by  electrons,  and  the  assumption 


made  that  the  ratio  of  this  scattering  in  the  superconducting  and  normai 
states  is  (T/T  )*  Tcf  equation  (19)1.  the  measurements  wouid  require  that 
H  .  (4  °K)  ^  at  ieas  t  iO'^  watt/crHdeg. ,  or  extra  resistance  wouid  be 

vislole  at  I  as  a  deviation  from  the  pure  behavior  observed.  Hence, 
tinea  X  .  -  has  been  observed,  at  lowest  temperatures,  to  be  limited  by 
boundary  Scatter ing,  this  gives  an  upper  limit  to  x  ^  ^  at  higher  temperatures 
which  Is  too  low  to  account  for  the  peak  in  the  Bufve  at  about  3  °K. 

Some  experimental  evidence^*  has  been  reported  which  seems  to  support  the 
second  interpretation.  Measurements  on  high  purity  tin  and  thallium,  Figs.  9 
and  10  show,  first,  a  behavior  of  the  electronic  thermal  conductivity  at  the 
transition  from  the  normal  to  the  superconducting  state  which  is  similar  to 
that  found  in  lead  and  mercury, and  second,  an  anisotropy  of 

'^o  "  -  0  °K’ 

which,  in  the  case  of  thallium  does  not  exceed  i0%,  but  for  tin,  cy  along 
the  rooil  *xis  is  i.4  times  greater  than  ry  along  the  TiiO]  axis,  see  Table  I. 

At  sufficiently  iow  temperatures,  the  electronic  thermal  conductivity 
in  the  normal  state  can  be  expressed 

T/h  -  p/L  +  T  W  (T)  .  (32) 

n  o  I 

Here  the  ffrst  term  is  determined  by  the  scattering  of  electrons  by  lattice 
defects  and  the  specimen  boundaries  (p  is  the  residual  electrical  resistivity 
and  L  is  the  Lorenz  number)  and  the  seBond  by  scattering  by  lattice  vibrations 
(phonons).  According  to  Makinson^’' 

TWj  (T)  .  n-T''  .  (33) 

There  is  also  a  change  in  the  quant i  ty 

-  dfT  W,(T)  l/dt''  ,  (34) 

which  on  the  simplest  theory''*'  should  be  a  constant.  The  change  in  T  W.  (T) 

(or  in  (vOt  see  Fig.  ii,  lies  outside  aii  possible  experimental  errors  and, 
can,  evidently  only  be  explained  by  a  lack  of  additivity  in  the  scattering  of 
electrons.^”  Table  I  shows  the  ratio,  at  T  «  T  ,  of  the  scattering  by 
lattice  imperfections  to  the  scattering  by  phonons,  which  is  determined  by 
the  value  of  p  /L  «  T^.  it  can  be  seen  from  the  table  that  for  the  purest 
specimens  measured,  the  conductivity  at  temperatures  near  T  is  limited  by 
scattering  of  electrons  by  phonons. 

The  data  of  Fig.  12  and  Table  I  show  that  for  a  relative  increase  in 
phonon  scattering  (a  reduction  In  o  /L  «  T ’’).  x  /x  decreases  near  T  . 

w  5  M  C 
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K  /h  is  then  close  for  specinens  of  thallium  and  tin"  with  the  same 
ratios  p  /L  oe  T  •’  (specimens  TX-3,  Sn-1  and  Sn-2) .  For  the  thallium  and 
tin,  therefore,  as  was  found  earlier"*"  in  the  experiments  on  mercury  and 
lead,  the  reduction  in  the  electronic  thermal  conductivity  at  the  transition 
from  the  normal  to  the  superconducting  state  is  sharper  when  the  electrons 
are  scattered  by  phonons  than  when  the  scattering  is  by  lattice  defects. 

A  qualitative  explanation  of  the  different  form  of  the  dependence  of 
on  T/T  for  scattering  by  imperfections  and  by  phonons  is:  In  the 
former  case, '"the  electron  mean  free-path  is  the  same  in  both  the  normal  and 
superconducting  states  of  the  metal,  in  the  latter  case  it  is  rather  the 
mean  collision  time,  t,  between  electrons  and  phonons  which  should  stay  con* 
stant,  assuming  an  unchanged  interaction  between  them,  a  result  of  the 
change  in  the  electron  energy  spectrum  in  the  superconducting  state,  this 
leads  to  a  reduction  in  the  electron  mean* free-path  and  therefore  in  the  value 
of  K  /k  .  A  calculation  of  the  heat  conductivity  for  the  case  r  ■  constant 
leadi  to  the  following  result 

~  t/T  ^  ~  **  2tp/expfAT"*  cosh  -1  +  1  ]drf  (35) 

o  o 

sech<l/2  ./  x’  +  (A/T)2ylx,  (36) 


where  ^  */  is  the  excitation  energy  in  the  superconductor, 

n  ■  (e  ' '  +  1)“^,  and  A  Is  the  superconducting  energy  gap,  assumed  to  be 
of  value  1.7  T  .  The  curve  R  of  Fig,  12  is  calculated  from  (35)  and  (36), 

It  can  be  seen^that  (35)  and  (36)  only  describe  the  variation  of  h  /k  with 
T/T  qualitatively.  This  may  possibly  be  due  to  the  electron-phonSn  Vnter- 
actfon  being  different  in  the  normal  and  superconducting  states  (by  electrons, 
is  meant  the  so-called  normal  electrons  of  the  superconductor).  Another  pos¬ 
sibility  Is  that  the  different  dependences  of  the  k  /h  on  T/T  for  thallium, 
tin,  mercury  and  lead  (Fig,  12)  are  produced  by  the*dl9ferencei  In  A/T  for 
these  metals.*^ 

The  Residual  Resistance  in  the  Superconducting  State 

in  the  case  of  tin,  indium,  tantalum,  thallium,  vanadium,  and  niobium, 
the  normal  : tate  electronic  thermal  conductivity,  for  those  specimens  which 
have  been  studied,  is  determined  at  T  and  below  by  imperfection  scattering, 
and  the  same  applies  to  the  various  afloys  and  to  impure  specimens  of  lead 
and  mercury.  These  specimens  can  be  used  to  test  (14)  for  k  ,  except 
where  k  ^  appreciable  and  complicates  the  picture.  As  TdAg  as  lattice 
_ ElLt5_ 

it 

There  appears  to  be  some  anisotropy  in  x^/h  vs  T/T  in  tin.  The  change 
of  K  /h  with  t/T  for  specimens  along  the  '’V lOlaxis^'fol  lows  a  steeper  law 
than*foP  specimens  along  the  fOOl]  axis.  This  is  similar  in  form  to  the 
anisotropy  In  ultrasonic  absorption  near  T 
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conduction  Is  unimportant,  the  ratio  k  /k  should  agree  with  (l4);  If 
K  /h  Is  larger  than  expected,  the  difference  Is  ascribed  though 

ll  l8  usually  not  possible  to  prove  that  this  Is  so,  ^  ’ 

K  has  been  measured  for  tin  down  to  about  T  /3.  Fig,  13  shows 

values®of  K  /h  the  specimens  are  numbered  in  Increasing  order  of  W  ,  Note 
that  the  cufveS  for  Sn2  and  Sn3  are  practically  coincident,  even  though  their 
values  of  W  differ  by  a  factor  of  about  2,  The  high  values  of  H Vx  are 
ascribed  to°k^  ,  Since  these  specimens  have  appreciable  lattice  tonauctlon 
even  In  the  nBnfial  state,  this  seems  a  plausible  Interpretation. 

Similar  results  were  obtained  for  indium.'^*  However,  in  all  cases 
it  is  found  that  the  observed  f- function  decreases  more  rapidly  with  ten>- 
perature  just  below  T  than  the  function  (14)  (semiempl rical  function).  It 
could  be  that  this  departure  from  the  conclusions  of  the  two>f1u1d  theory 
Just  below  T  Is  related  to  the  similar  departure  already  noted  for  the  g- 
functlon,  S?adek,  on  the  other  hand,  suggested  the  following  form  for 

’‘ejl,s^*‘e/,n‘ 

f  -  3(T/T  )V2  +  (TA ,  (37) 

c  c 


which  he  based  on  an  assumption  about  i  which  seems  no  more  artificial 
than  (13).*  * 


Substantially  similar  results  for  h  ^  have  been  obtained  for  tantalum®, 
tin.  Indium,  thallium,  tantalum  and  vanaBlum;'’*  the  case  of  niobium  was  com¬ 
plicated  by  frozen-ia  magnetic  flux. 


Not  much  data  arc  available  x  .  /h  ^  below  about  T  /2,  partly 
because  not  many  measurements  have  beeh’^de  sufficiently  fow  tenn 
peratures  and  partly  because  of  the  increasing  importance  of  k  .  as  the 
temperature  is  decreased.  ^ 


Measurements  of  n  below  about  I  °K  have  been  made.®^ » For  tin 
and  tantalum  indications  have  been  found^  that  h  /m  decreased  faster  with 
temperature  than  It  should  according  to  equations^( lS)  and  (l4).  This  was 
confirmed  by  the  extensive  work  on  tin^®,  wher#*  an  exponential  decrease  with 
temperature  was  found  in  k  for  two  pure  speci  changing  at  lower  tenw 

peratures  to  a  T®  dependence.  The  latter  variation  is  ascribed  to  lattice 
conduction,  while  at  higher  temperatures  k  .  is  important  (It  is  k  . 
which  varies  exponential  I  y),'‘®  A  speclmen^of^lead  which  had  prevlouSiy® 
bean  measured  at  helium  temperatures*”  was  measured  down  to  0.4  Below 

1  °K  It  was  found  that  k  «  T®  (presumably  lattice  conduction),  but  just  at 
the  upper  limit  of  the  measured  temperature  range.  Indications  were  found  of 
a  faster  variation,  which  Is  confirmed  if  one  joins  up  these  measurements 
with  those  at  higher  temperatures.  Probably  this  is  a  case  of  h  .  decreas¬ 
ing  exponentially  with  temperature.  Measurements  of  tin,  indium,  tBal 1 1  urn, 
aluminum,  tantalum,  and  niobium  from  0.4  °  to  1°  K  disclosed  an  exponential 
variation  of  in  the  case  of  thallium;  in  the  other  cases,  the  effect 
appeared  to  be  masked  by  lattice  conduction. 


14 


Presumably,  the  exponential  behavior  of  h  ,  ^  at  low  temperatures  is  re¬ 
lated  to  a  similar  variation  in  the  specific  hMt.  In  some  cases,  indications 
of  such  a  variation  of  the  specific  heat  of  superconductors  well  below  the 
transition  temperature  have  been  observed.  It  appears  likely,  therefore,  that 
the  Gorter-Casimlr  two- fluid  theory  breaks  down  at  very  low  temperatures. 

The  Lattice  Component  of  Thermal  Conductivity  in 
the  Superconducting  State 


In  superconductors  one  should  expect  heat  to  be  carried  by  phonons 
rather  than  by  electrons  at  low  temperatures,  where  the  concentration  of 
normal  electrons  must  become  vanishingly  small.  This  phonon  conduction 
will  be  enhanced  by  the  reduction  of  the  scattering  which  the  phonons  ex¬ 
perience  by  encounters  with  normal  electrons.  Evidence  for  this  effect  can 
be  obtained  either  by  observing  the  temperature  dependence  of  h  at  suf¬ 
ficiently  low  temperatures,  or  by  introducing  agents  into  the  specimens 
which  will  scatter  phonons  selectively.  Both  types  of  experiments  have  been 
carried  out,  and  have  shown  beyond  doubt  that  a  superconductive  metal  exhibits 
behavior  identical  with  that  of  a  dielectric  crystal,  as  far  as  thermal  con¬ 
ductivity  is  concerned,  when  sufficiently  near  to  absolute  zero.  Even  h  . 
is  too  small  to  be  observed  in  the  normal  state,  in  the  superconductive  sta?e 
decreases  very  rapidly  with  temperature,  so  that  at  very  low  temperatures 
H  ^  is  only  limited  by  the  scattering  of  phonons  by  static  imperfections  or 
bBuftia  ries.  Thus,  experiments  below  I  °K  usually  give  clear  Indications  of 
lattice  conduction. 

In  the  dimensional  equation. 


K  =  A  C  ;  V, 


(38) 


where  C  is  the  specific  heat  per  unit  volume  of  the  heat  carriers,  v  is  their 
velocity  and  }  is  their  mean-frec-path;  A  is  a  constant  which  is  usually  equal 
to  1/3,  V  is  temperature- independent  and  C  is  proportional  to  T^.’  Hence  h 
will  vary  with  T'  since  '  is  coni- tan t.  As  was  pointed  out  by  Casimir^”, 
scatter  of  phonons  only  occurs  at  the  geometric  surface  of  the  specimen  in  a 
perfect  dielectric  crystal  at  low  enough  temperatures.  This  means  that  " 
is  constant  and  is  determined  by  the  specimen  dimensions.  Thus  the  magnitude 
of  H  is  dependent  on  size  alone.  This  type  o*^  behavior  has  indeed  been  ob¬ 
served'®  in  dielectric  crystals,  the  heat  conductivity  being  proportional  to 
T"  and  of  the  predicted  magnitude. 

Clear  evidence  for  a  similar  behavior  in  superconductors  has  been  found 
in  measurements  of  thermal  conductivity  in  lead'-,  particularly  h  down  below 
1  K.  It  was  pointed  out'"  that  phonons  and  electrons  each  have  one  scatter¬ 
ing  mechanism  for  which  theory  and  experiment  agree,  and  one  for  which  there 
is  much  confusion.  The  heat  conduction  by  phonons  at  the  lowest  temperatures 
is  expected  to  be  limited  hy  boundary  scattering  and  to  vary  as  T’’ ,  This  is 
well  substantiated  experimentally;  Fig,  Il+  shows  results  for  lead  and  thal¬ 
lium.  Similar  results  have  been  obtained  for  tin,  indium,  and  niobium.  How¬ 
ever,  When  one  extrapolates  this  to  higher  temperatures  and  attempts  to  obtain 
that  part  of  the  phonon  conduc*"! ''i  ty  which  is  limited  by  electronic  scatter¬ 
ing,  confusion  arises  because  the  electronic  scatterers  are  themselves  car¬ 
riers  oF  heal.  As  a  result,  estimates  of  this  cart  of  the  phonon  conductivity 


are  very  unreliable  and  disagreement  with  theoretical  predictions  not  sur¬ 
prising;  similarly  with  electronic  conductivity  limited  by  phonon  scattering. 

Although  the  curve  of  Fig.  14  shows  a  definite  variation  of  with  T*, 
the  numerical  value  of  h  was  found  to  be  rather  smaller,  in  this  case,  than 
vK>uld  be  expected  on  the^basis  of  Casimir's  treatment  which,  for  a  cylindrical 
specimen,  predicts  /  »  d(d  is  the  specimen  diameter).  Similar  deviations  are 
found  In  other  cases.  At  the  lowest  temperatures,  k  has  been  found  to  vary 
as  in  lead^®,  tin  and  indium*’’,  and  tin.*’^*'’'  It^was  found  that  generally 
the  group  at  Oxford'"'*'^ » deduced  values  for  which  were  five  to  ten  times 
higher  than  values  expected  from  the  external  di^nsions  of  their  specimens, 
while  the  group  at  Cambridge*^*^^  found  values  of  Wg  in  rough  agreement  with 
the  external  dimensions,*  Two  possible  explanations  have  been  given'  for  the 
findings  of  the  Oxford  group:  (1)  their  specimens  contained  considerably  more 
grain  boundaries  and  (2)  is  reduced  by  frozen- in  magnetic  flux, 

A  study  of  the  effect  of  strain  on  k  shows  clearly  that  even  a  very 
small  degree  of  cold  work,  such  as  is  proiuced  by  a  slight  vibration  of  the 
cryostat,  can  raise  the  numerical  value  of  m  substantially.  This  difficulty 
was  not  fully  understood  in  most  of  the  reseirches  quoted  in  the  foregoing. 
More  careful  work*^  has  shown  that,  at  least  in  the  case  of  tin,  values  for  x 
can  be  obtained  which  are  sufficiently  close  to  the  value  predicted  by  Casimil 
to  make  scattering  by  the  specimen  boundaries  the  dominant  process. 

Even  so,  however,  the  tendency  of  h  to  fall  short  of  the  predicted  value 
deserves  attention,*  The  situation  can  oe  represented,  conveniently,  by  in¬ 
troducing  several  thermal  resistances  which  are  characterized  by  different 
mean-free-paths,  and  which  are  combined  addltively.  Thus,  denoting  the  mean- 
free-path  derived  from  observation  as  we  write* 

1//  -  (1/d)  +  (1/M  (39) 

O  D 

in  which  /.  is  a  mean-free-path  characteristic  of  an  additional  resistance  in 
the  bulk  material.  The  circumstance  tnat,  as  yet,  no  metal  specimen  has  been 
found  in  which  this  additional  internal  resistance  is  zero,  v/hereas  it  i s  of 
considerable  magnitude  in  most,  may  be  due  to  the  fact  that  the  metals  are 
much  more  plastic  than  the  dielectric  crystals  investigated.  This  is  a  field 
where  further  research  is  Indicated,  particularly  because  the  features  associ¬ 
ated  with  the  internal  scattering  appear  to  be  complex.  Although  the  tem¬ 
perature  dependence  of  h  is  cubic  in  many  cases,  both  higher  and  lower  powers 
of  T  have  been  observed.  Moreover,  it  has  been  found  that  even  at  the  lowest 
temperatures,  a  single  crystal  of  the  given  material  will  exhibit  a  thermal 
conductivity  proportional  to  T^,  even  though  a  polycrystalline  sample  may  show 
a  T*  dependence.  This  suggests  that  the  Internal  scattering  centers,  which 
are  most  probably  dislocations,  may  have  a  more  profound  effect  than  the 
crystal  boundaries. 

In  the  case  of  impure  tin  specimens**  h  varies  more  slowly  than  T*  except 
at  the  lower  end  of  the  temperature  range.  At  higher  temperatures  h  T'’. 
Goodman**  presumed  the  dominant  scattering  mechanism  In  the  T^  regioR"to  be 
scattering  by  the  free  electrons,  Klemens*  pointed  out  that  this  wa^  unlike¬ 
ly,  since  from  (19)  t  j  should  vary  as  T^,  not  as  T"’’,  However,  it  has  been 
shown*  that  dislocationS^can  play  an  important  part  in  determining  k  .  ,  even 
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in  normal  metals  where  Wp  is  not  quenched.  It  seems  quite  possible,  there¬ 
fore,  that  in  the  T”  region  x  .  is  limited  mainly  by  dislocations.  This 
interpretation  can  also  be  applied  to  a  tantalum  specimen*^  where  k  ^  ^ 
was  found  to  vary  as  T*  below  1  K.  * 

A  more  convincing  proof  of  the  phonon  nature  of  k  at  low  temperatures 
is  provided  by  work  in  which  selected  scattering  mechanisms  are  introduced 
into  the  specimen.'’  Observation  of  the  different  effects  of  changes  on 
and  K  allows  a  more  unambiguous  assessment  of  the  nature  of  the  energy 
transport  in  each  case.  Assuming  that  the  heat  is  entirely  carried  by  phonons 
in  the  superconductor,  at  low  enough  temperatures,  we  may  expect  the  following 
behavior  of  k  and  k  in  the  same  specimen: 

(l)  K  slioulS  be  insensitive  to  point  imperfections,  i.e.,  to  the  amount  or 
nature^of  the  impurity,  whereas  k  should  be  much  reduced  by  small  amounts  of 
impurity  and  should  be  sensitive  ?o  its  nature,  (2)  k  should  be  reduced  by 
large-scale  imperfections,  e.g,,  grain  boundaries,  such  as  are  introduced  by 
plastic  deformation  and  recrystallization,  v^liereas  x  should  not  be  greatly 
influenced  by  them.  Providing  that  the  specimen  is  pure  enough,  dislocations 
should  reduce  both  h  and  k  .  In  many  cases,  however,  k  is  so  small  already 
because  of  impuritiel,  that'^the  only  observable  effect  iS  the  reduction  in  k  • 
(3)  H  should  be  reduced  when  the  specimen  diameter  becomes  small  compared  ^ 
to  the^phonon  mean-free-path  in  the  bulk  material,  whereas  should  not. 

Experiments  designed  to  test  these  predictions  have  been  performed  on 
lead.^*  The  results  unequivocally  show  the  phonon  nature  of  the  heat  trans¬ 
port.  Fig.  15a  shows  and  k  for  two  single  crystals,  one  of  pure  lead  and 
one  of  an  alloy  containing  0.7?'  bismuth, “  The  conductivities  in  the  normal 
state  at  1  °K  differ  by  a  factor  of  100  and  (even  in  the  superconducting  state 
at  4  °K)  there  is  still  a  large  difference.  The  conductivity  at  this  tem¬ 
perature  evidently  is  still  mainly  electronic,  but  the  two  k  curves  merge 
below  2  °K.  Very  similar  numerical  values  of  H  were  obtalnid  at  1  °K  with 
specimens  of  lead  containing  the  same  amount  ortin  or  thallium.  On  the  other 
hand,  measurements  on  a  lead  specimen  containing  0,6%  thallium  {PbT,‘0,6%)  be¬ 
fore  and  after  severe  strain  due  to  bending  showed  h  reduced  to  1/6  its 
original  value,  whereas  h  was  completely  unchanged,  fig,  15b,®  Probably, 
the  dislocations,  introduced  into  the  sample  by  the  strain,  are  very  effective 
scatterers  of  phonons,  but  the  electronic  conductivity  is  unaffected  by  them 
because  of  the  tnallium  inpurity. 

The  effect  of  sample  size  on  k  is  more  difficult  to  demonstrate  clearly. 
Simply  to  compare  measurements  made  on  a  thick  rod  and  a  thin  wire  of  the  same 
material  can  not  be  conclusive  since  it  is  almost  impossible  to  avoid  strain¬ 
ing  the  latter.  Thus,  it  would  be  dilTicult  to  distinguish  with  certainty 
between  reduction  in  due  to  limitation  of  the  phonon  mean-free-path  result¬ 
ing  from  geometrical  boundaries  or  from  dislocations.  However,  an  experiment 
has  been  performed’"'  on  a  lead  foil  0.07  mm  thick  which  was  stabilized  me¬ 
chanically  by  being  rolled  into  a  scroll.  Since  the  phonon  mean-free-path  of 
the  material  was  of  the  order  of  0,5  mm,  a  size  effect  should  have  been  no¬ 
ticeable  in  the  foil.  Indeed  the  heat  conduction  of  the  foil  at  1  °K  was 
found  to  be  five  times  smaller  than  that  of  a  bulk  specimen  of  the  same  ma¬ 
terial,  whereas  k  is  essentially  the  same  in  both  cases.  Moreover,  the  tenw 
perature  dependence  of  k  for  the  scroll  approaches  T" ,  which  is  to  be  expect¬ 
ed  for  boundary  scattering. 


Metals,  such  as  the  transition  metals,  which  have  low  intrinsic  electrical 
conductivity  should  exhibit  the  most  pronounced  effect  of  enhanced  phonon  con¬ 
duction.  The  most  striking  example  of  this  kind  has,  indeed,  been  observed  in 
tantalum*'’  which  was  measured  down  to  0.2  °K,  The  result  of  these  investi¬ 
gations  Indicates  that  h  /h  Is  even  larger  than  one  in  the  neighborhood  of 
1  If  the  data  are  plotted  in  a  form  which  permits  comparison  with  the 

semlempl rical  function  (f-function) ,  see  Fig.  16,  it  can  be  seen  that  the 
latter  represents  the  data  remarkably  well  down  to  about  0.4  T^.  Below  this 
temperature,  there  is  an  enormous  rise  in  the  phonon  conduction  wh>ch  at  0.2  T^ 
assumes  a  value  about  a  hundred  times  larger  than  the  electronic  contribution. 
Data  on  niobium.  Included  In  Fig.  16,  and  on  vanadium  present  a  similar  pat¬ 
tern.  The  results  for  tantalum  have  been  used  to  separate  the  phonon  contri¬ 
bution  K  from  the  electronic  part  h  .  This  analysis  shows  that  at  a 

temperatBr4®T  /4,  the  phonon  conductivity’ is  Si^ll  increasing  with  decreasing 
temperature  even  though  the  ratio  h  .  /m  ,  ^  approaches  10*',  indicating  the 
strong  scatter  of  phonons  by  the  normST  conduction  electrons.  Comparing  the 
numbers  clearly  indicates  that  the  maximum  in  k  will  be  less  pronounced  or 
will  disappear  for  metals  with  better  intrinsic  electronic  conductivity. 

Aluminum  is  a  typical  example  of  the  latter  behavior.  Here  the  data 

agree  well  with  the  semiempirical  function  (f-function)  over  the  entire  range. 
Including  the  lowest  temperatures.  Fig,  17  shows  the  observations”-^  of  h 
plotted  for  pure  aluminum  and  an  alloy  with  1,7  atomic  %  copper.  The  purl 
metal  follows  the  B.C.S.  function  with  an  energy  gap  slightly  smaller  than 
that  predicted  by  the  theory.  At  temperatures  belov;  T  /3,  however,  the  alloy 
shows  a  deviation  to  higher  values,  indicative  of  a  phSnon  component. 

The  thermal  conductivity  has  been  measured  on  two  indiunv- lead  alloy 
samples  containing  4,05  and  7.31  atomic  %  lead,  respectively,  down  to  0,4 
The  results  were  interpreted  in  the  light  of  the  Bardeen,  Rickayzen,  and  Te- 
wordt  (BRT)^'’’  theory.  This  theory  assumes  that  the  dominant  electron  scatter¬ 
ing  Is  due  to  impurities  and  the  dominant  phonon  scattering  is  caused  by  e- 
lectrons.  The  normal-state  data  were  fitted  with 


=  AT  +  BT^  (40) 

where  the  first  term  represents  k  -  and  the  second,  k  .  ,  This  determined 

the  two  parameters  A  and  B,  It  wis’?ound  that  the  phonBn’Sontribution  was 
very  small  in  both  the  normal  and  the  superconducting  cases.  Similarly  the 
superconducting  data  were  expressed  by 


H  ■  ATR  +  H  . 
s  ex  ph,s 


fUi) 


Here  R^^  Is  the  BRT  ratio  of  electronic  thermal  conductivity  '  the  two 
states,  and  k  ^  the  lattice  component,  included  the  effects  of  boundary 
and  point-defect  Scattering  as  well  as  the  electronic  scattering  considered 
by  BRT.  Hence,  the  additive  resistance  approximation 


'^*‘ph,s  "  '^**BRT 

was  used.  Here 


(42) 
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(43) 


''•ORT 


=  BT' 


is  tiie  BRT  value  assuning  phonon  scattering  hy  electrors  alone,  Kgp  is  the 
phonon  thernai  conductivity  iinited  by  boundary  and  poini-defect  scattering 
according  to  Siac!  ;ipp  involves  tv#o  additional  adjustable  paraneters,  the 
mean- free- pa  til  dve  to  boundary  scattering  L  and  a  temperature  \'/hicli  is  a 
measure  of  the  point-defect  scattering,  in  evaluating  L,  the  velocity  of 
sound  was  taken  as  1,0x10"  cn/sec.  The  accepted  value  of  the  energy  gap  for 
indium,  ^^^(Oy/Jk  =  0.7  was  used. 

As  shown  in  Figs,  1"  and  i9  the  agreement  betiiKien  theory  and  experiment 
is  quite  good.  The  fit  to  the  superconducting  data  obtained  in  this  manner 
is  particuiarly  good  at  the  iower  temperatures.  In  this  region  boundary  and 
point-defect  scattering  predominate,  Tiie  slight  disagreement  near  1  n  is 
attribi:ted  to  the  limitations  o'  the  additive  resistance  approximation,''"  The 
small  phonon  contribution  ( liic  linearity  of  the  normal-state  data  rules  out  the 
strong  possibility  o'  a  larger  plionon  term)  probably  accounts  for  the  poor  fit 
at  higher  temperatures,  v/here  the  BRT  result  predominates.  However,  it  is 
quite  apparent  that  k  is  largely  lattice  conductivity  over  almost  the  entire 
range  shown.  The  nah^num  or  plateau  region  near  I  “it  for  the  superconducting 
case  is  caused  by  the  transition  from  the  scattering  of  the  phonons  by  elec¬ 
trons,  '.1^,  ‘o  a  combination  of  boundary  and  poInt-dfii'ccL  scattering, 

One  errnot  obtain  a  value  for  V.'  from  the  very  lov/  temperature  ob¬ 
servations  (e<'<cept  an  upper  limit,  wnf^h  is  probably  very  much  larger  than 
Wg  ),  because  o''  tne  importance  o'  phonon-phonon  interactions,  which  are  not 
influenced  by  the  superconducting  behavior.  It  is  only  possible  to  observe 
V/_  O'er  a  limited  temperature  range  below  T,.,  It  has  been  determined  in 
this  way  for  Sn96-Hg  lead-bismuth  alloys''  which  had  been  previously 
measured’’",  and  “^or  i ndium- thal  1  lum  al  loys,  " 


A  ronf  1  i ct i I’ll  pictwire  c-  the  ratio  h  •=  V/  /Wg  5  results  from  these 
observations,  Huir*  suggested  h  .v  (T  /T)'  ,  0TS?n'''’'’h  •».  (T  /T)'',  and  Sladek'’", 
whose  measurements  seen  most  sui table^for  the  evaluation  o^  h,  did  not  obtain 
a  simple  po’-/er  'cvi,  nor  the  sane  curve  of  n  ’'ersus  T/T  for  al  i  his  samples, 
but  rather  a  series  of  cur'/es  for  i,  nil  in  the  vicinity  of  h  =  (T  /J)',  but 
too  high  juS"  below  t  and  ^endin.g  to  becoi.ic  too  low  at  io'«*r  temperatures. 


The  i  ntnrpre !  at  ion  contains  uncor  t  ai  .-.t  i  os ,  hence  ,  could  easily 
have  conponeiits  than  (.his  is  certainly  so  at  ?ov;esl  te.mperat.iros 

and  nay  '  e  so  even  If  ,  <-,r  T"  ) ,  wuicli  v;ould  tend  to  decrease  h.  The  sep¬ 

aration  of  in*^o  H  .  ^'and  ;t  ,  in'ol'es  the  assumption  that  the  ratio  f  = 

H  /h  is  indcpernefif  o''  alTbt'  connosiMon.  Nevertheless  there  are  pro- 
bS.'^^y  ?4al  discrepa  ncies  'ro  •  I'  =  (T  /T)  ,  part i c.;iarly  just  below  T  ,  just  as 
them  are  discrepancies  in  >•.,  ^  ; us  1 '  :  eio'.-.'  T^, 


!  \ios  le.i ;  lone  I  cbo'  o 
breaks  lo'.'n  in  .iosr  cases 


C-ai  :he  '-'onctior,  ( ser  I emp i  r i ca i  formula) 
li  './liicu  sca'ii.ring  o'  electrons  in  the  normal 


sia.c  is  1,0  caused  impur'l-ies  tw  .  •  p.’ionons.  In  'act,  tl.e  observed 
v'alues  o'  r<^/u  'or  a  metal- li'ce  lead,  which  has  a  lov;  characteristic  tem- 
pprature  and  a  relatively  ■in,  tronslion  temperature  (7.2  °K) ,  do  not  follow 
the  '-''TClon  for  any  terporoture  region.  1'  can  be  seen  'rom  Fig.  1 
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however,  that  exhibits  a  maximum  at  0,5  T  which  is  similar  to  those  shown 
by  tantalum  and  niobium  at  much  lower  reduced  temperatures.  The  question 
therefore  arises  whether  the  maximum  in  »  of  lead  is  different  in  nature 
from  the  enhanced  phonon  conduction  found^in  the  h  of  the  transition  metals. 
Fortunately,  the  sensitivity  of  the  phonon  conduction  to  dislocations  allows 
this  problem  to  be  decided  experimental ’y,**"*  When  the  pure  single  crystal 
of  lead  of  Fig,  1  was  strained  at  helium  temperature,  it  was  observed  that  y. 
was  Indeed  reduced.  However,  tliis  reduction  occurred  only  at  temperatures 
well  below  the  maximum  in  >t  ,  i.e,,  in  the  reduced  temperature  region  found 
in  tantalum  and  niobium. the  maximum  in  lead  was  quite  unaffected,  whereas 
It  was  drastically  reduced  in  the  transition  metals  by  similar  treatment. 

A  clearer  picture  of  the  behavior  in  the  case  of  lead  can  be  obtained 
from  the  curves  In  Fig.  20.^  This  shows  that  at  temperatures  above  ~  0.4  T 
the  experimental  values  are  lower  than  those  given  by  the  f-function 
(semlempi rical  function),  whereas  they  are  higher  below  that  temperature. 
lntroduc!ng  dislocations  into  the  specimen  by  strain  does  not  affect  the 
higher  temperature  region  of  k  at  all.  However,  k  now  follows  the  semi- 
emplrrcal  function  remarkably  wll  below  ~  0,4  T  ,  ^Thus,  it  seems  that  the 
failure  to  obey  the  f-function  arises  from  two  quite  different  reasons.  At 
high  temperatures,  where  the  predominant  process  is  the  scatter  of  electrons 
by  phonons,  the  theoretical  understanding  Is  not  yet  sufficient,  whereas 
phonon  conduction  becomes  predominant  below  0.4  T  .  Once  phonon  conduction 
Is  drastically  reduced  by  scattering  on  dislocations,  the  semiempi rical 
formula  holds  quite  well.  Hence  the  maximum  in  h  in  lead  is  entirely 
electronic  in  nature  and  is  clearly  connected  witS  the  maximum  in  y 

n 

KIcmens*  points  out  that  the  interpretation''’’ of  the  thermal  conductivity 
of  lead-bismuth  alloys  contains  difficulties.  Fig,  21  shows  plots  of  k  and 
K  versus  T  for  these  alloys.  It  is  easily  seen  that  if  the  increase  o?  k 
fSr  the  alloys  0,2%  Bi  and  0.5%  Bi  over  the  values  of  k  for  alloys  of  low^ 
bismuth  content  Is  to  be  explained  in  terms  of  enhanced^lattice  conduction, 
then  H  .  for  these  two  alloys  is  higher  than  h  .  for  the  alloys  0.1%  Bi 
and  0,Bz%*Bi,  and  possibly  even  higher  than  h  *’”for  pure  lead.  Of  course, 
it  is  possible  to  explain  this  by  assuming  soRe'tmperfect ions  to  be  present 
in  the  more  dilute  alloys  and  not  in  the  more  concentrated  alloys,  though 
this  disagrees  with  the  usual  observations. 

Thermal  Conductivity  in  the  Intermediate  State 

A  superconducting  pure  element,  in  the  shape  of  a  long  cylinder,  under¬ 
goes  a  sharp  transition  from  the  superconductive  to  the  normal  state  upon  the 
application  of  a  longitudinal  magnetic  field.  In  other  cases  the  transition 
Is  gradual;  Increasing  the  magnetic  field  causes  a  gradual  Increase  of  flux 
In  the  specimen,  until  all  the  material  is  in  the  normal  state.  Upon  removing 
the  field,  the  material  does  not  return  to  the  original  superconductive  state, 
but  some  magnetic  flux  remains  frozen  in. 

The  Intermediate  state  of  the  material  is  not  homogeneous,  rather  it  con¬ 
sists  of  a  mixture  of  normal  and  superconducting  regions  in  the  material,  the 
former  having  high  flux  density  (above  the  critical  field)  and  the  latter,  zero 
flux.  Since  the  lines  of  flu;:  are  continuous,  the  structure  of  the  intermediate 
State  Is  dominantly  one  of  filaments  or  layers,  alternately  normal  and  super¬ 
conducting,  lying  in  the  direction  of  the  field. 
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A  number  of  measurements  have  been  made  of  the  thermal  conductivity  of 
superconductors  in  the  intermediate  state,  with  the  specimens  as  long  cylin¬ 
ders,  With  longitudinal  fields  there  will  usually  not  be  a  marked  mixing  of 
•^ne  two  states,  except  In  the  case  of  alloys,  when  the  normal  state  inclusions 
'■'ill  be  mainly  filaments  running  the  length  of  the  specimen.  With  transverse 
“^ields,  however,  the  specimen  will  readily  break  up  into  a  mixture  of  two 
states,  and  the  normal  state  inclusions  would  then  be,  predominantly,  layers 
■'erpendlcular  to  the  cylinder  axis,  and  thus  to  the  direction  of  heat  flo’/.’, 
f->e  thickness  of  the  Individual  regions  nay  be  of  the  order  of  0,01  cm. 

With  a  cylindrical  specimen  in  a  longitudinal  field,  with  normal  and 
superconducting  filaments  along  the  direction  of  heat  flow,  the  overall  thermal 
conductivity  would  be  expected  to  be  given  by  the  average,' 

K  =  V  K  +  ( 1  -  ■■:  )  n  i^) 

n  n  ns 

•■•here  v  is  the  fraction  of  normal  material  and  can  be  deduced  from  flux  meas- 
urement§.  Similarly,  for  transverse  fields,  the  thermal  resistance  averages 
are  given  by 

U  =  V  W  +  (1  -  V  )W  ,  (45) 

n  n  n  s 

As  a  consequence  of  either  (Ab)  or  (45),  the  thermal  conductivity  in  the 
intermediate  state  (either  with  subcritical  fields  or  with  frozen-ln  (trapped) 
flux)  should  be  intermediate  between  the  normal  conductivity  h  and  h  as 
measured  in  the  virgin  superconducting  state.  Within  the  boun9s  of  t^is  re¬ 
striction,  hysteresis  effects  are  possible. 

Such  behavior  was  indeed  observed  in  many  early  measurements  in  transverse 
fields'’*®®  and  in  longitudinal  fields,'  Later  measurements  showed  variations 
of  K  with  magnetic  field  strength,  however,  which  could  not  be  reconciled  with 
either  (44)  or  (45),  In  the  case  of  some  lead-bi smuth  el  loys  end  of  niobium, 
it  has  been  found'"’  that  the  thermal  resistance  passed  through  a  maximum  value, 
on  applying  a  field,  which  was  higher  than  the  resistance  in  either  the  normal 
or  the  superconductive  state.  Fig,  22,  Also,  on  removal  of  the  field,  x  did 
not  return  to  the  original  value  x  ,  but  to  a  value  again  lower  than  either  x 
or  K  (a  hysteresis  effect).  These  anomalies  seem  more  likely  to  occur  the 
loweP  the  temperature  and  the  larger  the  impurity  content,  but  this  Is  not  a 
general  rule.  This  effect  has  been  ascribed’*'’  to  a  heat  How  mechanism  which 
occurs  in  the  superconductive  state  In  addition  to  electronic  conduction;  In 
the  Intermediate  state  this  mechanism  would  be  inhibited,  A  two-fluid  circu¬ 
lation  was  Suggested  for  this  additional  mechanism,'*''  A  more  likely  explan¬ 
ation  seems  to  be  lattice  conduction,''' 

Further  instances  of  anomalies  in  transverse  fields  were  found  later  for 
pure  lead'’'*’ »'*■’ ,  pure  tin  and  indium  (but  not  for  tin  containing  0,134%  Bl)"®, 
and  for  mercury,"'’*'"’ 

The  anomalously  high  thermal  resistance  in  the  intermediate  state  has 
been  attributed  to  various  scattering  mechanisms  acting  at  the  boundaries  of 
the  filaments.  The  situation  has  been  greatly  clarified  by  a  systematic 
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investtgation  of  :ii':  alloy  sfirie*;  Snir..'"  The  induced  magnetic  moment  and  the 
dependence  of  the  electrical  and  thermal  resistance  on  the  magnetic  field 
strength  were  measured  on  these  specimens.  A  typical  set  of  results  for  an 
almost  perfect  single  crystal  containing  2,8%  indium  is  given  in  Fig.  23  for 
two  temperatures.  Whereas  the  fraction  of  frozen- in  flux  Is  largely  tem¬ 
perature  independent,  the  maximum  in  the  thermal  resistance  is  pronounced 
at  2,03  and  practically  nonexistent  at  2,6$  °K,  From  this  it  can  be  de¬ 
duced  that  the  structure  of  the  intermediate  state  is  much  the  same  at  dif¬ 
ferent  temperatures  and  that  the  relevant  scattering  mechanisms  may  be  dif¬ 
ferent,  More  detailed  analysis  does  indeed  suggest  that  the  filament  boundar¬ 
ies  scatter  both  electrons  and  phonons  but  that  the  temperature  dependence  is 
different  for  the  two  mechanisms. 

Theory  of  the  Thermal  Conductivity  of  Superconducting 
Alloys  with  Paramagnetic  impurities 


The  prediction  of  "gapless"  superconductivity  in  paramagnetic  alloys'"' 
and  the  confirmation  of  this  prediction*'*'  are  significant  recent  develop¬ 
ments,’*  Although  the  theory  is  based  op  an  approximeitc  treatment  of  a  simple 
model  and  more  detailed  experiments  are  needeo,  the  prediction  is  unambiguous 
and  the  confirmation  convincing.  The  assumption  is  that  the  static  magnetic 
impurities  are  randomly  distributed  and  that  their  spins  are  uncorrelated. 

The  theory  makes  it  clear  that  the  key  feature  of  the  superconducting  state 
is  the  condensation  phenomenon.  A  gap  In  the  single  particle  excitation 
spectrum  is  evidently  not  a  necessary  requirement  for  either  the  infinite 
conductivity  or  the  perfect  diamagnetism  of  the  condensed  state.  The  experi¬ 
mental  investigation  of  this  phenomenon  of  superconductivity  without  an  energy 
gap  promises  to  Improve  our  understanding  both  of  superconductivity  and  of  the 
effects  of  magnetic  impurities  in  metals  in  general. 

Based  upon  this  theory,  the  electronic  thermal  conductivity  h  of  a  weak¬ 
ly  coupled,  isotropic  superconductor,  doped  with  a  small  concentration  of 
paramagnetic  impurities,  is  calculated,  starting  from  a  Kubo  formula,  by  con¬ 
sidering  the  electron- impuri ty  interaction  in  the  ladder  approximation. 

A  considerable  simplification  of  the  final  expression  occurs  if  it  is 
assumed  that  the  total  single  particle  lifetime  Is  much  smaller  than  the  ex¬ 
change  lifetime  t  .  With  this  assumption,  an  expression  is  obtained  for  H  /h 
of  the  form  ^ 

CO 

Hs/Hn  ■  (3/2tt^)3®  r  dw  u/*  sech”(i^'=(w)h((ju/A,»)  (46) 

o 

or  »  ,  (47) 


A  is  the  average  order  parameter,  and 


3  =  (1/kT) 


(48) 


Is  is  the  Boltzmann  constant,  and  T  absolute  temperature.  The  units  are  chosen 
so  that-fi  =  1,  For  non-magnetic  irpurilies,  the  usual  result  is  found  that 


h(u)/A  <  1)  =  0  ('i?) 

and 

h(u/A  >  1)  ■  1  ,  (50) 

the  average  order  parameter  a(T,t  )  being  the  energy  gap  in  this  case.  For  a 

paramagnetic  alloy,  the  lower  limit  of  integration  is  the  physical  energy  gap 

uu  (a,(v)  and  not  the  average  order  parameter  A(T,t  ),  Moreover,  the  function 

h?<u)  Increases  smoothly  toward  unity  for  «,  >  lu  , 

o 

The  ratio  k  /k  is  evaluated  numerically  as  a  function  of  the  reduced 
temperature  (T/T  )  »  t  for  different  impurity  concentrations,  including  non¬ 
magnetic  impurities.  The  results  are  shown  in  Fig,  2h, 

Abrikosov  and  Gor'kov  have  shown  that  the  energy  gap  function  u:  (T)  Is 
quite  different  than  the  Gor'kov  order  parameter  \(T)  in  such  alloys,®  It  is 
found  that  v  /h  is  less  than  unity  in  the  "gapless"  region  At  <  1,  This 
theory  predilts^that  even  in  the  gapless  region  the  thermal  conductivity  In 
the  superconducting  state  is  lower  than  that  in  the  normal  state,  because  al¬ 
though  the  energy  spectrum  has  no  gap,  it  is  stili  distorted.  In  addition, 
the  onset  of  gapiessness  does  not  lead  to  an  abrupt  change  in  the  thermal  con¬ 
ductivity,  Long  before  the  gap  actually  vanishes,  the  BCS  singularity  in  the 
density  of  states  is  smoothed  out  by  the  impurity  scattering.  Finally,  it  is 
found  that  h^/k  ps  a  function  of  t,  has  a  characteristic  concentration  de¬ 
pendence,  For  ?  ~  0,9,  H  decreases  for  small  n.  while  for  t  '^0,75,  It 
Increases  with  hj(nj  is  the  paramagnetic  impurity  concentration)* 

For  non-magnetic  Impurities  In  systems  with  weak  electron- phonon  inter¬ 
actions  (to  which  this  theory  is  restricted),  a  simple  caiciiiatio.i  of  the 
thermal  conductivity,  using  a  Boltzmann  equation,  is  possible.  To  justify 
the  elaborate  formalism  used,  it  is  pointed  out  that  in  the  gapiess  region 
(always  close  to  T  ),  it  is  not  possible  to  associate  a  narrow  band  of  en¬ 
ergies  with  a  state  of  momentum  near  the  Fermi  momentum,  i,e.,  the  quasi- 
particie  approximation  breaks  down.  In  addition,  the  effects  of  the  para¬ 
magnetic  impurities  in  renormalizing  the  energy  spectrum  are  crucial.  Both 
these  effects  would  cause  difficulties  in  conventional  transport  theory,  but 
they  are  easily  tal'.en  into  account  in  the  Kubo  formulation. 


The  Thermal  Conductivity  as  a  Defect  Detector 
in  Superconductors 


It  has  been  demonstrated  aljove  that  rhonon  conduction  predominates  in  k 
at  low  temperatures,  that  it  is  sensitive  to  strain,  and  independent  of  * 
point  imperfections.  Clearly  tnis  can  provide  o  dis'^inction  between  these 
two  different  types  of  lattice  defects.  Investigations  of  this  nature  have 
ttie  added  advantage  that,  in  the  same  specimen,  the  scattering  of  phonons 
and,  by  simply  applying  a  magnetic  field,  that  of  electrons  can  be  studied 
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without  even  warr,!,.j  ‘.ne  sample  o  roon  ter.perature.  In  most  cases  the  mag¬ 
netic  fields  required  to  destroy  percoiiductivi ty  are  less  than  1000  oersteds. 
The  accompanying  magneto- resistive  effects  in  the  normal  state  are  small. 

Those  v/hic:-.  exist  can  usually  be  taken  into  account  by  running  comparison  ex¬ 
periments  above  T^,'* 

Systematic  experiments  lave  been  performed  on  the  effect  of  strain  and 
Impuri ty,‘'^»‘^‘'‘  The  clearest  results  were  obtained  by  Rowell‘S*,  who  subjected 
pure  lead  and  a  lead  alloy  to  controlled  bending  at  helium  temperatures  and 
measured  the  thermal  conductivity  in  the  normal  and  superconducting  states 
before  and  after  Introducing  the  strain.  Another  feature  was  the  study  of 
annealing  effects  at  different  temperatures  between  that  of  liquid  helium  and 
room  temperature.  Order  of  magnitude  agreement  was  obtained  between  the 
density  of  dislocations  derived  from  the  measurements  of  heat  conductivity  and 
those  predicted  on  the  basis  of  the  strain  Introduced.  More  spectacular  re¬ 
sults  than  were  obtained  on  lead  are  those  on  a  niobium  rod  which  was  origin¬ 
ally  single  crystal  and  which  was  subsequently  stretched  in  steps  until  it 
ruptured,  h  showed  no  effect  except  at  the  highest  strain.  Even  after 
fracture  x  changed  only  by  a  few  per  cent.  On  the  other  hand,  the  effect 
on  H-  is  fSr-reaching.  In  its  undisturbed  condition,  the  niobium  specimen 
Showed  a  very  pronounced  maximum  in  x  at  temperatures  below  0.4  T  .  With 
successive  stretching,  this  maximum  was  largely  removed.  Point  imperfections, 
which  scatter  electrons,  will  reduce  both  x  and  h  at  high  reduced  temper¬ 
atures,  where  the  heat  transport  is  still  b9  elections.  The  maximum  in  k  at 
low  temperatures  because  of  phonon  conduction  is  unaffected.  However,  since 
x_  Is  drastically  reduced,  x  may  now  exceed  x  in  this  temperature  region, 

Ac  still  lower  temperatures, *th Is  relation  is  again  reversed  and  k  /h  becomes 
smaller  than  unity.  Extended  lattice  defects,  such  as  dislocation!,  ?ntro- 
duced  by  strain,  appear  to  have  little  or  no  effect  on  the  electronic  part  of 
the  conduction  mechanism.  Therefore,  x  ,  as  well  as  the  high  temperature  part 
of  X  t  will  not  change  materially.  In  ?he  low  temperature  region,  on  the 
othe?  hand,  where  heat  is  carried  by  phonons,  h  is  decreased.  The  phonons 
are  scattered  by  the  dislocations;  with  Increasing  strain  the  ratio  h  /k 
becomes  progressively  smaller,^  ® 

This  method  can  also  be  applied  to  the  study  of  nuclear  radiation  damage. 
The  thermal  conductivity  of  a  single  crystal  rod  of  niobium  was  measured  first 
in  the  undamaged  state  then  after  neutron  irradiation  at  room  temperature. 

Both  K  and  x  are  reduced  by  irradiation.  Fig.  25.  If  it  is  assumed  that 
both  iRterstlclals  and  vacancies  have  been  produced  by  irradiation,  it  is  be¬ 
lieved®  that,  although  the  former  may  have  migrated  at  the  temperature  of  ir¬ 
radiation,  the  temperature  was  never  high  enough  to  cause  migration  of  va¬ 
cancies  in  niobium.  Hence  the  decrease  in  x  is  ascribed  to  the  vacancies 
produced  directly  by  the  radiation  and  assume  a  more-compi icated  process  to  be 
responsible  for  the  decrease  in  x  .  A  calculation  based  on  the  condensation 
of  interstitials  due  to  irradiation  yields  an  increase  of  about  10'’  disloca- 
t lons/cm®  in  niobium.  Using"^" 

0  -  (h^v^b^/28  k')N  (51) 

for  a  random  array  of  dislocations,  where  h  is  Planck's  constant,  v  the 
velocity  of  phonons,  v  is  Gruneisen's  constant,  b  the  magnitude  of  the 
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durt^er'b  vector,  k  8o!tzmann  constant,  and  N  Is  the  density  of  dislocations 
per  unit  area,  to  determine  the  change  In  k  with  the  number  of  dlsj^ocatlon 
lines,  analysis  ot  the  observed  change  In  h  yields  3x10^  llnes/cm",  which  Is 
In  surprisingly  good  agreement.  ^ 

•  ,  Conclusions  ■' 

The  thermal  conductivity  of  superconductors  is  discussed  from  the  view¬ 
point  of  the  two-fluid  model  and  the  Bardeen-Cooper-Schrieffer  theory.  The 
significance  of  the  ratio  y  /h  is  discussed  and  its  dependence  on  T/T  Is 
shown  to  be  a  universal  function  independent  of  the  particular  element^meas- 
ured.  The  electronic  thermal  conductivity  of  superconductors  with  strong 
electron- phonon  coupling  differs  markedly  from  typically  weak  superconductors 
like  tin  or  indium.  A  theory  is  discussed  to  expiain  this  phenomenon.  The 
thermal  conductivity  in  the  superconductive  state  Is  discussed  from  the  stand¬ 
point  of:  (a)  the  ideal  resistance,  (b)  the  residual  resistance,  (c)  the 
lattice  component,  and  (d)  in  the  Intermediate  state.  The  existence  of 
"gapless"  superconductivity  is  shown  in  the  theory  of  the  thermal  conductiv  Ity 
of  superconducting  alloys  with  paramagnetic  Impurities.  Finally,  the  use  of 
thermal  conductivity  to  detect  defects  in  superconductors  Is  discussed,  t 
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rhe  thermal  conductivity  of  (a)  tin  and  (b)  lead  in  the 
normal  and  superconducting  states.” 
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Fig.  19.  TheriMl  conductivity  vorsus  temporaturo  for  ln(7«3lVb) 


Fig,  20,  Th«  •ff^ct  of  ttroln  on  /m  vorsus  T/T  for  lood.® 
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Flq.  24.  Th*  r«tlo  of  the  electronic  thermel 
conductivity  In  the  superconducting  and  normal 
states^  ’‘*^*'n  ''*  *  *  T/T^  for  various  paramag¬ 
netic  Impurity  concentrations.  Is  the  tran¬ 
sition  temperatura  for  the  rcisvant  Impurity 
concentration,  the  latter  being  expressed  In 

terms  of  n^,^,  the  concentration  required  to 
complataly  destroy  superconductivity.  The 
curve  for  nonmagnetic  Impurities  alo’«e  Is  de¬ 
noted  by  Oj  •  O.OOn^^.  The  n|  ■  O.TOn^^  curve 
Is  almost  Identical  to  the  0.85  curve  for 
t  ~  0.7  and  Is  not  shown  explicitly  In  this 
region.  These  are  based  on  equation  (46).^* 
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by  etching  specimen  Tl-7. 


DOCUMENT  CONTROL  DATA  •  R  &  D 

'Sgruf^fy  cims'titication  ot  thie,  body  ol  abstract  mc'emr**  ti  I'totatton  mutt  hf  gnfgrgd  ivhan  the  oyaral*  report  is  e laaaltiad) 


\  OHIOiNATtNO  ACTli/tTv  t  CorpofAfg  author  j  *^«.BePORTSEC^JPlTy  CLASSIFICATION- 

^  !  Unclassified 

U.S.  Army  Electronics  Command 

Fort  Monmcjth,  New  Jersey 


3  REPORT  title 


THERMAL  CONDUCTIVITY  OF  SUPERCONOUCTORS 


4.  DESCRIPTIVE  NOTES  (Type  ol  report  and  inclusive  dalea) 

Technical  Report 


9.  Au  Tho  RiS>  ('Firsf  n«m«.  middle  initial,  la  at  name' 


Irving  N,  Greenberg 


REPORT  DATE 

December  19^>9 


7a.  total  '.o  or  PACES  7^.  NO  OP  PEPS 

54  76 


OISTrirUTIOS'  SIVTC'CHT  (D 

Tlili  ^otwr  gnl  hos  ItCTi  public 

filiQSP  eng*  tote;  lit  diiftibwtion  l>  vnlimlfrd. 

n  SuRPL £MEN T ART  NOTES 

li  SPONSORING  MILtTAPY  ACTIVITY 

institute  for  Exploratory  Research 
AMSEL-XL-G,  U.S.  Army  Electronics  Command 
Fort  Monmouth,  New  Jersey  07703 

This  report  Is  one  of  a  series  which  Is  intended  to  inform  the  reader  on  the  pre¬ 
sent  status  of  the  thermal  conductivity  In  solid  materials.  The  first  of  the  series 
was  USAELROL  Technical  Report  2361,  dated  May  1963,  which  gave  a  general  overall  picture 
of  the  field  but  specifically  discussed  the  conductivity  of  dielectric  solids  (Insul¬ 
ators).  This  was  followed  by  Technical  Report  ECOM-2799.  dated  February  1967,  which 
dealt  with  concepts  mainly  applicable  to  metals,  such  as:  a.  The  electronic  component 
of  thermal  conductivity  and  Its  relation  to  the  free  electron  theory,  b.  The  telation- 
ship  between  electrical  and  thermal  conductivity,  as  exemplified  by  the  Wiedraann-Franz- 
Lorenz  law,  and  the  limitations  of  the  law,  c.  Electron- lattice  wave  scattering,  and 
the  restrictions  on  the  theory.  Next  In  the  series  was  Research  and  Development  Tech¬ 
nical  Report  ECOM-2932,  dated  February  1968,  wherein  considerations  mainly  applicable 
to  semiconductors  and  semimetals  were  discussed,  such  as: 

a.  The  Influence  of  carrier  concentration  and  degeneracy  on  the  thermal  conductivity  o( 
both  extrinsic  and  intrinsic  semiconductors,  and  on  the  WIedmann-Franz-Lerenz  law. 

b.  The  electronic  component  of  the  thermal  conductivity,  In  terms  of  both  holes  and 
electrons,  for  an  extrinsic  semiconductor,  c.  The  transport  properties  of  an  Intrin¬ 
sic  semiconductor,  d.  The  thermal  resistance  in  semiconductors  arising  from  free 
electrons  (holes). 

After  this,  Research  and  Development  Technical  Report  ECOM-3026,  dated  September  1968, 
was  written.  This  treated  the  lattice  component  of  thermal  conductivity  In  metals. 
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13.  Abstract  (Cent) 


alloys,  and  semiconductors  and  the  Influences  affecting  it.  Phonon  scattering 
processes  ware  described  by  a  relaxation  tl>ne. 

2  2 

Since  aii  intrinsic  resistivities  are  proportional  to  C*  a.C.  ,  and  since 

J  J 

the  coefficients  C.  cannot  be  obtained  fro^  first  principles,  the  lattice 

component  was  expressed  in  terms  of  the  intrinsic  electronic  component  of 
thermal  conductivity.  The  residual  thermal  resistivl ty  was  also  considered. 
The  intluenca  of  the  ide^l  or  intrinsic  electronic  thermal  resistivity  and 
that  of  the  residual  thermal  resistivity  on  the  total  thermal  conductivity 
tMS  discussed,  A  rcther  detailed  discussion  was  given  of  the  methods  of 
separating  the  measured  thermal  conductivity  into  <he  electronic  and  lattice 
components.  The  intrinsic  lattice  component  of  thermal  conductivity  and  in¬ 
fluences  affecting  it  were  considered  rather  completely.  A  fairly  exhaustive 
treatment  of  the  effect  of  crystal  imperfections  on  the  i.nttice  thermal  con¬ 
ductivity  wes  given. 

In  the  present  report  we  discuss  the  thermal  conductivity  of  super¬ 
conductors  from  the  viewpoint  of  theory  set  forth  in  the  previous  reports 
plus  considerations  from  the  theory  of  superconductivity. 
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